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Abstract 

One of the continuing challenges in cosmology has been to determine the large-scale space-time 
metric from observations with a minimum of assumptions - without, for instance, assuming that 
the universe is almost Friedmann-Lemaitre- Robertson- Walker (FLRW). If we are lucky enough 
this would be a way of demonstrating that our universe is FLRW, instead of presupposing it or 
simply showing that the observations are consistent with FLRW. Showing how to do this within 
the more general spherically symmetric, inhomogeneous space-time framework takes us a long way 
towards fulfilling this goal. In recent work researchers have shown how this can be done both in the 
traditional Lemaitre-Tolman-Bondi (LTB) 3 + 1 coordinate framework, and in the observational 
coordinate (OC) framework, in which the radial coordinate y is null (light-like) and measured down 
the past light cone of the observer. 

In this paper we investigate the stability of solutions, and the use of data in the OC field 
equations including their time evolution - i.e. our procedure is not restricted to our past light 
cone - and compare both approaches with respect to the singularity problem at the maximum of 
the angular-diameter distance, the stability of solutions, and the use of data in the field equations. 
We also compare the two approaches with regard to determining the cosmological constant 
A. This allows a more detailed account and assessment of the OC integration procedure, and 
enables a comparison of the relative advantages of the two equivalent solution frameworks. Both 
formulations and integration procedures should, in principle, lead to the same results. However, 
as we show in this paper, the OC procedure manifests certain advantages, particularly in the 
avoidance of coordinate singularities at the maximum of the angular-diameter distance, and in 
the stability of the solutions obtained. This particular feature is what allows us to do the best 
fitting of the data to smooth data functions and the possibility of constructing analytic solutions 
to the field equations. Smoothed data functions enable us to include properties that data must 
have within the model. 

PACS numbers: 98.80.-k, 98.80.Es, 98.80.Jk, 95.36,+x 
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I. INTRODUCTION 



Over the last 15 years, there has been great progress made in cosmology. The concor- 
dance model - a flat FLRW universe with Q m « 0.27 and Qa ~ 0.73, where these fi's 
represent the matter content and dark-energy content of the universe, respectively - has 
received continuing strong support from an explosion of detailed observational data together 
with its careful analysis vis-a-vis theoretical considerations. This model seems to provide 
an increasingly reliable account of all the observations relevant to the structure and history 
of our universe. 



However, there are lingering uncertainties and unresolved controversies - particularly 
with regard to dark-energy (is it simply the cosmological constant A, something dynamical, 
evolving, or an alias for the influence of intermediate and large-scale inhomogeneities?) and 
with regard to the overall adequacy of perturbed FLRW models in a very inhomogeneous 
universe, at least on small and intermediate length scales. It would be reassuring to actually 
demonstrate from observations that the metric of the universe on scales larger than a certain 
definite value - or averaged over volumes larger than a certain definite value - is almost 
FLRW, instead of simply assuming that, and then finding that the observations are more 
or less consistent with almost-FLRW. One of the difficulties, of course, is that they might 
also be consistent with one or more cosmologies which are far from FLRW. At present 
we have no assurance that an almost-FLRW universe is the unique best-fit model to the 
observations, even on the largest scales. We do not even know reliably what is the smallest 
scale on which the universe is almost-FLRW. And it is obvious that it certainly is not on 
small and intermediate scales. 



There are certainly strong indications that our universe is almost-FLRW, as already 
emphasized. But in order for these to constitute definitive evidence more is needed - in par- 
ticular either observational evidence that the cosmological principle (Copernican principle) 
is correct, that we are not located in a large scale under-dense bubble, for instance, or that 
observational data constraining the models is almost-FLRW. For, a careful critical treatment 
of these issues, see jl|. As they emphasize, any such conclusions demand norKDerturbative 



investigations and confrontation with the data. Though some recent studies [2 



44[ argue that 



data on spectral distortions due to scattering of CMB photons off free electrons between 
z — 1 and z — 10 significantly restrict the kind of void which may surround us, more precise 
and extensive measurements together with non-perturbative modeling beyond FLRW is 
essential, as Clarkson and Maartens [l] stress. 



In order to begin to provide a broader framework for eventually banishing these uncer- 
tainties - and possibly providing confirmation of the large-scale almost-FLRW character of 
our universe - there has been considerable recent theoretical work oriented towards showing 
how we can begin from an assumed cosmological ansatz considerably more general than 
FLRW, and determine its metric from observations in a rigorous way. This is possible in 
principle [5], though meshing the observational and the theoretical has proven very tricky, 
and accumulating enough precision data of th e rig ht sorts is still uncertain. In particular, 



a number of workers, including ourselves, |6l-ll9| have begun with a general spherically 
symmetric cosmology and then allowed the cosmologically relevant observational data to 
determine the solutions to the field equations. This has been done using the usual 3 + 
1 LTB coordinates and showing in detail how the observations can be used to determine 
the space-time metric by numerical integration of the field equations. Others, including 
ourselves, have shown, alternatively, how to integrate the general spherically symmetric 
field equations in observational coordinates (OC), in which the observer's past light cone 
plays a prominent role. In fact, the radial coordinate in this framework is a null parameter 
measuring "distance" down the generators of the observer's past light cone. 



In this paper, we shall discuss how the discrete observational data, with its gaps and 
uncertainties, can be handled - either by keeping it discrete and binning it by redshift or 
by smoothing it into "data functions"; whether there are any difficulties with coordinate 
singularities at the maximum of the angular-diameter (observer area) distance, which is 
one of the key components of the observational data; and whether or not there are stability 
problems with the equations themselves. Will small uncertainties in the data lead to a "blow 
up" of the solutions we are seeking? Careful analysis of these aspects of the OC integration 



procedure will allow us to compare and contrast it with the 3+1 LTB approach [13Hl5|. 
and to assess the advantages and disadvantages of each. At the same time, we shall compare 
how the data can be used - particularly the maximum of the angular-diameter distance - 



in the two approaches to determine the value of the cosmological constant A for the more 
general LTB - not just for FLRW - models. 



In the next section we define observational coordinates, write the general spherically 
symmetric metric using them and present the central conditions for the metric variables. 
Section II I II presents the basic observational parameters we shall be using and several key 
relationships among the metric variables. Section IIVI presents the full set of field equations 
for the spherically symmetric case, with dust and with A / 0. In Sec. |V] we present a general 
integration scheme for all inhomogeneous spherically symmetric LTB models. Section |VT] is 
dedicated to the comparison between the the OC formulation and integration and that of 
the usual (3 + 1) approach. In Sec. IVIII we examine the stability of the solutions to the null 
Raychaudhuri equation. Since that is the first step in solving the field equations with data 
on our past light cone, stability of its solutions is necessary for insuring the stability of the 
solutions to the other field equations. This complements the discussion of stability of the 
OC integration in the last paragraphs of Sec. IVI1 Sec. IVIIII discusses the implementation of 
an algorithm for using the spherically symmetric OC field equations with data to determine 
the metric. Section [IX] is devoted to our conclusions. 

The new results in this paper are found in Secs JVIl IVI!] and IVIIII They include an 
explicit and detailed comparison of the OC and the 3 + 1 LTB approaches to integrating 
the spherically symmetric field equations, a discussion of advantages and disadvantages 
of solving the equations numerically in contrast to first fitting the data with smooth "data 
functions" , a brief explanation of how the 3 + 1 approach can use the maximum of the angular 
diameter distance to determine the cosmological constant A, a study of the stability of the 
OC formulation of the field equations, particularly the null Raychaudhuri equation, and a 
detailed example (in Sec. IVIIip of carrying out the full integration of the OC equations, 
using FLRW data functions. Here no stability problems are encountered. Sections [TT] - 
IVl provide a succinct, but detailed synthesis of the OC approach, the OC field equations, 
and their integration with data functions on the light cone, and then off it to obtain the 
time-dependence of the space-time metric functions. 
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II. THE SPHERICALLY SYMMETRIC METRIC IN OBSERVATIONAL COOR- 
DINATES 



In observational coordinates {w, y, 6, <p} the Spherically Symmetric metric takes the gen- 
eral form: 

ds 2 = -A(w, yfdw 2 + 2A(w, y)B{w, y)dwdy + C(w, y) 2 dtt 2 , (1) 

where w is the time coordinate defined such that w =constant specifies a single past light 
cone along the observer's world line C. Our past light cone, here and now is given by w = wq. 
y is the null radial coordinate measured down the null geodesies generating each past light 
cone centered on C, which is given by y — 0. y increases as one moves down the past light 
cones from C and is chosen to be co-moving with the matter flow, that is u a d a y = 0, u a being 
the fluid 4-velocity satisfying u a = A^d^ and u a u a = —1. 9 and are simply the latitude 
and longitude of observation - i.e. spherical coordinates based on a parallely propagated 
orthonormal tetrad along C. 

Besides the freedom to choose y co-moving, there are other coordinate freedoms which 
preserve the form of the metric Eq. ([1]). In particular we can choose w to be any time 
parameter we like along C - that is, at y = 0. This is usually done by choosing A(w,0). 
Furthermore, we can choose y to be any null distance parameter on an initial past light cone 
- usually on w = wq. Since y is co- moving, that choice is then dragged off by the fluid flow 
onto other light cones. We shall use this freedom to set 



For furt 
Stoeger [U 



A(w ,y) = B(w ,y). 
rer details about observational coordinates see Ellis et al 



12| 



(2) 

[s| and Araiijo and 



It is important to specify the central conditions for the metric variables A(w,y), B(w,y) 
and C(w,y) in Eq. ([!]) - that is, their proper behavior as they approach y = 0. These are: 



as y -> : A(w, y) ->• A(w, 0) ^ 0, 
B(w,y)^B(w,0)^0, 
C(w,y)^B(w,0)y = 0, 
C y (w,y)^B(w,0). 



(3) 



These conditions insure that C, our world line, is regular - so that all functions on it our 
bounded, and that the spheres (w, y = const) go smoothly to C as y — > 0. They also 
insure that the null surfaces w = const are past light cones of observers on C (See Ref. j^] , 
especially section 3.2, p. 326, and Appendix A for details). 

There are several other relationships which are important for later use. In particular we 
have two fundamental four-vectors in the problem, the fluid four-velocity u a and the null 
vector k a , which points down the generators of past light cones. These are given in terms 
of the metric variables as 

From the normalization condition for the fluid four-velocity, we can immediately see that 
it can be given (in covariant vector form) as the gradient of the proper time r along the 
matter world lines: u a = —r, a . It is also given by ([1]) and @ as 

u a = g a bU b = -Aw ja + By :a . (5) 

Comparing these two forms implies 

dr = Adw — Bdy <^ A = r w , B = —r y , (6) 

which shows that the surfaces of simultaneity for the observer are given in observational 
coordinates by Adw = Bdy. The integrability condition of Eq. ([6]) is simply then 

A' + B = 0. (7) 

where a "dot" indicates d/dw and a "prime" indicates d/dy. This turns out precisely to be 
the momentum conservation equation, which is a key equation in the system and essential 
to finding a solution. 

III. THE BASIC OBSERVATIONAL QUANTITIES 

The basic observable quantities on C, which are consistent with the central conditions, 
are the following: 
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(i) Redshift. The observed redshift z at time wo on C for a comoving source a null radial 
distance y down C~(po) is given by 

A(w a ,y) 

(ii) Redshift-drift (time-drift of the cosmological redshift). In principle we can measure 
the change in the cosmological redshift of distant galaxies from one time of observation to 
another (as a function of redshift). This is the redshift drift z(z). As a function of y, using 
Eq. (jSJ), it can be written 

■I \ (~\ _i \M( w o,0) i(w Q ,2/)l 



As many others have shown, and we have confirmed within the OC framework, [ 12|, |20 



221 ]; see also Appendix |A] ] z contains evolution- and galaxy-count-independent informa- 
tion about the mass-energy density. With reliable measurements of z and of observer-area 
(angular-diameter) distances, we can determine the mass-energy density h(wq,z). This 
would enable us to avoid having to use galaxy number counts (see below) to determine 
that. It turns out, besides, that i as a measurement on our past light cone can be used 
to determine the important relationship z = z(y), which enables us to translate key data 
as functions of z into functions of y [see Appendix |A] . Without z data, as we shall see 
below, z = z(y) can be determined only by solving the null Raychaudhuri equation [Eq. 
(122"]) below], which needs, instead, reliable galaxy- number-count and average galaxy mass 
data, which is very difficult to obtain. 

Though we do not yet have the technological capability to measure z, the prospects 



for dev elop ing it over the next 20 years are good. Loeb 23| and more recently Pasquini 



et al. 



241 ] have studied in depth the requirements and the promising ways of doing so 



using the hydrogen and metal absorption lines in Ly-a forests. Pasquini, et al. [24J using 
simulations of the spectra along with analysis of the expected noise, as well as the projected 
precision and stability of the planned European Extremely Large Telescope (E-ELT) and 
its CODEX ultra-stable spectrograph, conclude that a measurement of redshift drift with 
significant precision would be possible with these instruments over a ten-year time-line [ see 
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also Dunsby et al. 



25| 



(iii) Observer area distance (Angular diameter distance). The observer area distance, 
often written as r , measured at time Wq on C for a source at a null radial distance y is 
simply given by 

r = C(w ,y), (10) 

provided the central condition (j3J), determining the relation between C(w,y) and B(w,y) 
for small values of y, holds. This quantity is also measurable as the luminosity distance di 



because of the reciprocity theorem of Etherington 



26| (see also Ellis |23|), 



d L = (l + z) 2 C(w ,y). (11) 

(iv) The maximum of observer area distance. Generally speaking, C(wp, y) reaches a 
maximum C m .g X for a redshift z max between about 1.6 and 2.0 [Hellaby 
and Tivon 



; see also Ellis 



29| and Araiijo and Stoeger [30J ]. At C max , of course, we have 
dC(w ,z) dC(w ,y) 



dz dy 

further conditioned by 

d 2 C(w , z 



0, (12) 



dr> < °' (13) 
These C max and z max data provide additional independent information about the cosmology. 

Without C max and z max we cannot constrain the value of A with our other data. 

(v) Galaxy number counts. The number of galaxies counted by a central observer out to 

a null radial distance y is given by 

N(y) = 4vr f fi(w , y)m- 1 B{w Q , y)C{w , yfdy, (14) 
Jo 

where /i is the mass-energy density and m is the average galaxy mass. Then the total 
mass-energy density can be written as 

dz 1 

H(w ,y) = mn(w ,y) = M (z) — — -, (15) 

dy B{w ,y) 



where n(wo,y) is the number density of sources at (wo,y), and where 

mil dN 

J dfl rl dz 



Here dVt is the solid angle over which sources are counted, and J is the completeness of the 
galaxy count, that is, the fraction of sources in the volume that are counted is J. The effects 
of dark matter in biasing the galactic distribution may be incorporated via m and/or J . 
In order to effectively use number counts to constrain our cosmology, we shall also need an 
adequate model of galaxy evolution. We shall not discuss this important issue in this paper. 
But, fundamentally, it would give us an expression for m = m(z) in Eqs. (fl5|) and (fl6|) 
above. 

In line with our discussion of redshift drift earlier in this section, Araujo and Stoeger 
[12I have recently shown within the OC framework that the mass-energy density of the 
Universe can be fully determined in terms of the cosmological redshift, the time drift of the 
cosmological redshift, and the observer area distance [see Appendix |A]. Using the redshift- 
drift as a replacement for mass-energy density element in the minimally required data set 
needed to construct an LTB model for the Universe means that the number counts data is 
no longer needed. Moreover, as pointed out by Liske et al. {31 ] and Linder {^J one of the 
main features of the redshift drift measurement is its model independency. 

In identifying, choosing and discussing these basic observable cosmological quantities 
or probes in the general spherically symmetric OC context, we are not intending to ex- 
clude the potential important contribution of other probes that are being studied and show 
promise -including weak gravitational lensing, baryonic acoustic oscillations, various cosmic 
microwave background radiation (CBR) measurements, and CBR-related measurements of 
the integrated Sachs-Wolfe( ISW) and the Sunyaev-Zel'dovich (SZ) effects. These essen- 
tially provide observational avenues, like those we have discussed here, to determining the 
mass-energy density of the universe, angular-diameter distances, etc. We have not included 
discussion of them here, because they are not as straight-forwardly related to cosmological 
parameters in the general OC context, and often rely, at least in many applications so far, 
on assuming that the universe is FLRW - or perturbed FLRW - on the largest scales. 



IV. THE SPHERICALLY SYMMETRIC FIELD EQUATIONS IN OBSERVA- 
TIONAL COORDINATES 



Using the fluid-ray tetrad formulation of the Einstein's equations developed by Maartens 
33( and Stoeger et al [34( , one obtains the spherically symmetric field equations in observa- 
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tional coordinates with A ^ (see Stoeger et al [6] for a detailed derivation). Besides the 
momentum conservation Eq. (J7J), they are as follows: 

A set of two very simple fluid-ray tetrad time-derivative equations: 



(B C\ 

u, = -3^ + ^Y (18) 



C\ 6 

where \x m again is the relativistic mass-energy density of the dust, including dark matter, 
and 

1 C C 1 / C \ 2 
^y) = -^ + ACBC + 2{Bc) ' (19) 

is a quantity closely related to f-i mo (y) = Mm(w>o,y) [ see Eqj28l below]. 

Eqs. (IT7|) and ffTSl can be quickly integrated to give: 

B(w ,y) C 2 {w ,y) 

Hm{w,y) = Vm (y) r 7^7 — r; 20 

B(u;,y) C 2 (w,y) 

u(w, y) = (u (y) + —j - - = -— + — — + - ^— j , (21) 

where co (y) = u(w ,y) and the last equality in (I2T]) follows from the definition of u given 
above. In deriving and solving these equations, and those below, we have used the typical 
A equation of state, pa = —fJ-A, where pa and /ia = are the pressure and the energy 
density due to the cosmological constant. Both uq and /io are specified by data on our past 
light cone, as we shall show, /ia will eventually be determined from the measurement of 
C m ax and z max . 

The fluid-ray tetrad radial equations are: 



c 



C \ A B 



(co (y) + f)C 3 (w ,y) 



-\Hmo B(w ,y) C 2 (w ,y) 



c 



+ 



c 



QL — M.QL 

c b c 



(22) 
(23) 
(24) 
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The remaining "independent" time-derivative equations given by the fluid-ray tetrad formu- 
lation are: 



^ = S4 + (^ + ^M 2 ; (25) 



C CA V 2 

B B A , r, 1 , r, , 

B = BA ~ 2U A 2^ A - (26) 



From Eq. ( )23|) we see that there is a naturally defined "potential" (see Stoeger et al jsj) 
depending only on the radial null coordinate y - since the left-hand-side depends only on y, 
the right-hand-side can only depend on y: 



y W 2 J + I> (27) 



Thus, from Eq. (|23l itself 



^o(y) = -if - 7^7 7 / Hmoiy) B(w ,y) C 2 (w ,y) F(y) dy. (28) 

6 2C- S {w ,y) J 

All these spherically symmetric equations in observational coordinates have been indepen- 



dently obtained by Hellaby and Alfedeel 16] directly from the metric, without relying on 
the fluid-ray tetrad equations. 

We now introduce a simple but very important observational relationship originally ob- 
tained by Hellaby [28] in the 3+1 framework [for a detailed derivation of this result in 



observational coordinates see Araujo and Stoeger [111 

6M max + fi A C 3 max - 3C max = 0. (29) 
where, M max is the maximum of the quantity M(y) given by 

M(y) = -L F mN'(y)F(y)dy = -L F M'(y)F(y)dy. (30) 

87T JO 87T JO 

where, M(y) = mN(y) is the total mass summed over the whole sky by a central observer 
out to a null radial distance y. 

Eq. ( 12 9 p has to be considered a fundamental relation in Observational Cosmology, since 
it enables, from C(wo, z max ) and z max measurements, the determination of the unknown 
constant jjlj^. 
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Finally, from Eq. (I14p we have 



mN'(y) = 47ifi mo (y)B(w ,y)C 2 (w ,y). (31) 
Hence, Eqs. f )30|) and fl3Tj) give 

M(y) = \ [\m (y)B(w ,y)C 2 (w ,y)F(y)dy. (32) 
2 Jo 

This clearly shows the relationship between M(y) and the mass-energy density. 

V. THE GENERAL SOLUTION - TIME EVOLUTION OFF OUR LIGHT CONE 

In this section we outline the general integration procedure that is applicable to all inho- 
mogeneous spherically symmetric universe models - that is the only constraint. [See Araiijo 
and Stoeger [ll| for a detailed description of this procedure.] We do not know whether the 
Universe is homogeneous or not. But the data give us redshifts z, observer area distances 
(angular diameter distances) r (z), the angular-diameter distance maximum C max (w , z) at 
Zmaxi and galaxy number counts, which together with an overall average mass per galaxy 
can in principle give us Mq(z). Eventually with the advent of very large telescopes and 
very high precision spectrographs, as we have already indicated, we may also have z(z) 



data, giving us the time drift of cosmological redshifts |23_|, |24| . It is important to specify 
CmaxiwQ, z) at z max , because, as we have already emphasized, without them, we do not have 
enough information to determine all the parameters of the space-time in the A ^ case. 
Although we can determine C(wq, z) with good precision (by obtaining luminosity distances 
d>L and employing the reciprocity theorem, Eq. ( TTTT) ) out to relatively high redshifts, at 
present we do not yet have reliable data deep enough to determine C max and z ma x- But this 
is now within the realm of possibility, with precise space-telescope distance measurements 
for supernovae la, and the potential for the linear-size of ultra-compact (miliarcsecond) 
radio sources as a standard rod, enabling reliable C(wq, z) out to redshifts z ~ 4 (such data 



already exists, 
substantiated) 



3ut the status of these sources as standard measuring rods must be better 



35|. 



In pursuing the general integration with these data, we use the framework and the inter- 
mediate results we have presented in Sec. [TV] Obviously, one of the key steps we must take 
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now is the determination of the potential F(y), given by Eq. (I2~T|) . This means we need to 
determine C'(wo,y) and C(wo,y), which we now write as C and Co, respectively. We also 
need A(w ,y). We remember, too, that on w = w we have chosen B(w ,y) = A(w ,y), 
which we have the freedom to do. 



Clearly, C' can be determined from the r (z) = C(wo,z) data, through fitting, along 
with the solution of the null Raychaudhuri equation, Eq. (122]) . which has the first integral 

l^* 1 ^) 2 ^! 1 ^!* 1 ^^^ 5 }"' (33 > 

to obtain z = z(y) upon inversion (Stoeger et al p). As already mentioned, this enables 
us to write all of our data as functions of y instead of as functions of the redshift z. In 
Sec. IVll and particularly in Sec JVII] we shall discuss this important step in the integration 
procedure, focusing on the stability of its solutions. We determine Co by solving Eq. (124"]) 
for w = w . Its solution, taking into account the appropriate boundary condition (central 
condition) is given by: 



My) 



l 



My) 



2A' (y)C (y)C (y) 



{C'MY ~ My)C*{y)»A dy. (34) 



2C (y) Jo V UV " A (y) 
This procedure enables us to determine F(y), which obviously also depends on n^. Our 
next step is to insert this result along with N'(y) - the spatial variation of the galaxy number 
counts - and the average mass per galaxy m into Eq. ( 130]) to determine the mass function 
M{y). Next we evaluate the mass function M{y) at y max and plug the result into Eq. (|29|) 
which becomes an algebraic equation for fi\. With this determination of fi\, we know C (y) 
completely, and can now determine F(y) from Eq. (1271) . Furthermore, we observe from Eq. 
(TiT?]) that the quantity Uo(y) is also completely determined at this stage. From here on, we 



can now follow the solution off w = wq 
It is shown in Araujo and Stoeger 
following equation for A(w ,y): 



or all w. 



ll| that from Eqs. ([7]) and (1261) one obtains the 



My) = My)< 



2 ^o (y) + ^ mo (y) 



A 2 (y) dy + d 



(35) 



where we have written Ao(y) and Ao(y) for A(wo, y) and A(wo, y) respectively. Since we can 
set A(w,0) = 1 [see Appendix [B] , it is obvious that C x = A (0)/A (0) = 0. 
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Hence, we have shown that the data on our past light cone determines Ao(y). Since we 
know A Q (y) from the data, Eqs. ( 125]) and (1261) evaluated on our past light cone become 
algebraic equations for Co(y) and B (y), respectively 



- C Q {y)A Q (y) 
°ol2/J — TT~\ h 

A'MMv) 



Ao 2 (y)C (y) 



B (y) 



Mv) 



2uj o (y) + ^moiy) 



A 2 (y)B (y) 



(36) 
(37) 



where in the later we have used Eq. ([7j). 

Araujo and Stoeger have shown that the next step in this procedure is to obtain an 
equation for A (w ,y) through differentiation of Eq. (126]) with respect to w, and use of Eq. 
(J7|). That leads to: 



A {y) idw 



2u + -Hm) BA- 



}dy\ - A (y)C 2 , (38) 



where, for the same reason as above, C 2 = A (0)/A (0) = 0. 

It is important to note that all quantities on the right-hand side of the above equation are 
obtainable either directly from the data or from the algorithmic steps in the procedure we 



are outlining here [see Araujo and Stoeger 



111 ] for a detailed description of the algorithm]. 



Therefore, we have shown that we can obtain Ao(y) from the data. It is clear now that 
repetition of this procedure will give us all time derivatives of A, B and C on our past light 
cone, which means that A(w,y) B(w,y) and C(w,y) are completely determined by data on 
our past light cone, and calculable as Taylor series. 

We see from the above procedure that each step begins by finding the successive time 
derivatives of the metric function A(w,y) on our past light cone, d™A(wo, y), which are 
completely determined given our choice of the time coordinate w as measuring proper time 
along our world line C. 
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VI. COMPARISON OF THE OBSERVATIONAL-COSMOLOGY (OC) AND OR- 
THOGONAL (3 + 1) (MH) INTEGRATIONS 



The formulation of general spherically symmetric cosmologies in OC coordinates we have 
mst presented (and discussed more fully elsewhere 

HQ) 

is equivalent to that which Bondi 
3(1, Bonnor 37], and more recently Mustapha, Hellaby and Ellis 13], Lu and Hellaby [jjj 
and McClure and Hellaby [15J have developed in orthogonal (3 + 1) coordinates, where the 
radial coordinate is spacelike, instead of null, as it is in our OC treatment. Integration of 
the equations with data in this earlier alternative formulation (which we shall from now on 
refer to as the Mustapha-Hellaby (MH) formulation) of LTB models has been developed 
and discussed in detail by Mustapha, Hellaby and Ellis (hereafter MHE) [3], and in much 
more detail by Lu and Hellaby 14] and McClure and Hellaby fisl ] . Here we briefly sketch 
key aspects of that integration and compare it with the one we have just discussed, empha- 
sizing the common features and the differences. This comprehensive comparison of the two 
approaches is one of the primary aims of this paper, and is new to the literature. 



The MH formulation appears to be quite different to the equivalent OC formulation. 
This stems principally from the fact that the MH form of the metric is expressed in terms, 
as we have already mentioned, of a spacelike radial coordinate, which is orthogonal to the 
time coordinate, whereas in OC the radial coordinate is null, and measures distance down 
successive past-light cones of the observer. Thus, MH sports a foliation of spacelike surfaces 
normal to the temporal axis, while OC is characterized by a foliation of past light cones, 
which are, of course, null surfaces, along the observer's world line - a past light cone for each 
instant of time. As mentioned in the OC literature, the principal reason for this choice is to 
reflect the fact that all the observational data we receive is really arrayed on the observer's 
past light cone. Furthermore, all the cosmological fluid particles - astronomically speaking, 
the galaxies, like the observer's own galaxy - move along world lines which intercept the 
past light cones at different locations along it - that is, at different redshifts and at different 
angular positions on the plane of the sky. So, the OC formulation is more natural. More 
importantly, transforming from OC to MH coordinates, or vice-versa, presents difficulties. 
There is in general no exact solution of the null geodesic equation for relating a spacelike 
radial coordinate to a null radial coordinate. There are also a number of other very conve- 
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nient features of the OC formulation (see below), despite some conceptual challenges. 



A less obvious difference between the two formulations and integration procedures is that 
in the MH framework the emphasis is on solving for the three arbitrary functions in the 
system of field equations and exploring how the evolution of the space-time is determined 
by them (see [3] and references therein). In the OC approach the emphasis is on using the 
observational data directly to determine the metric functions, whose precise relationship to 
the data is highlighted and worked out before carrying out the integration. Although MHE 



recommend fitting the relevant data with smooth data functions, [13j in the most recent 
detailed implementations of the MH case [lsj the approach is discrete, reflecting the 
discreteness of the data and the need for numerical integration of the field equations. 



In the OC approach, in contrast, though there is no explict avoidance of a numerical, 
discrete approach, we and our collaborators have envisioned having analytic functional forms 
which the discrete data obeys, with the possibility of carrying out an analytic integration of 
the field equations, at least in some simple cases. This is driven by the realization that the 
observer-area-distance-redshift and galaxy-number-count (or, equivalently, the mass-energy- 
density- redshift) data for an exact FLRW universe, whether with or without a cosmological 
constant A, do take theoretically derivable forms with definite characteristics - and can be 
analytically integrated, as we show here and elsewhere [see [7, 9, 10 1 and references therein]. 
In fact, as has been also shown, a space-time is FLRW if and only if that data can be exactly 
fit by one or other of those functional forms [see [7|, [lOj and Sec. IVIIII below]. In practice 
this approach, of course, necessitates fitting suitable smooth analytic functions to discrete 
data, after correcting them, and keeping track of the deviations of the actual corrected data 
from the data functions. Another reason for doing this is to illustrate more clearly the 
actual integration algorithm - the sequence in which the steps in the algorithm must be 
implemented - and to provide a check on possible later numerical integration procedures. 
Then if some or all of the steps in a given case need to be numerically integrated, that can 
be done against that informed context. 



Even with these strong differences, however, there are clear similarities, as there must be. 
Among these, and one of the most prominent, is the need to solve the null Raychaudhuri 
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equation, Eq. (|22|) . and more immediately its first integral Eq. ( 133|) - in order to obtain 
z = z(y) or z = z(r) for the OC and the MH approaches, respectively [r is the spacelike 
radial coordinate in the MH formulation]. This key relationship is needed, because all the 
data is given as functions of the redshift z, and must be rewritten as functions of either 
y or r, respectively. The solution scheme for this equation is almost the same in both 
formulations [compare for instance Eq. (34) of Stoeger, et al. Q with Eq. (51) of MHE [3] 
- notice also that in neither version is there a singularity in the solution of this equation, 
at the maximum of the observer-area distance (see below)]. However, later work by Hellaby 
and coworkers somewhat obscures this similarity. Interestingly enough, too, in both cases 
(generalizing the MH treatment to the A ^ case) this equation is independent of the 
cosmological constant A. This facilitates the integration. 



With determination of z — z(y), or in the MH case z = z(r), the data functions for 
Cq(z) = r (z), N(z) or M (z), and in the future z(z), can, as already indicated, be written 
as functions of y or of r on our past light cone in each approach. This is another obvious 
feature common to the both formulations. However, in the OC case knowledge of z = z(y) 
immediately enables us to write down the g 00 component of the metric, A (y) = A(w , y), on 
our past light cone w = Wq by Eq. ([8]). The freedom we have in choosing y also enables us 
to write B(wo,y) by Eq. ([2]), and the data then gives us the other metric variable on wq as 
well, Co(y) = C(wo,y) = r (z(y)). Thus we quickly have the solution of the field equations 
on our past light cone, w = Wq. The next and more difficult part of the integration is to 
move this solution off wq into the past. 



Accomplishing this requires, as seen in the two previous sections of this paper, deter- 
mining the time-derivatives of the metric variables on wq, as well as F(y) and M(y). Once 
we have z = z(y) from Eq. (133|) . and have written the data as functions of y, we solve Eq. 
flM|) for Co(y). This enables us to calculate F(y), and then M(y), via Eqs. f )27|) and fl30l) . 
respectively. We already know B (y) from Eq. ([7]) B (y) = — A' (y) . As we have also 
just seen, Ao(y) can be determined via Eq. ( 1351) . and higher time derivatives of A, B and 
C are found through Eqs. ( 1381) . (1371) . ( 1361) . and their higher-order generalizations. These 
enable us to construct a Taylor-series solution in the time w for all the metric variables. In 
carrying this out, there are no inherent calculational difficulties or singularities. C'(y) = 
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at y = Umax, but this quantity never appears in a denominator - and so by itself never 
induces a singularity or an instability at y ma x- Furthermore, the actual equations we are 
employing to find the solutions are at this stage linear. Once we know M(y), of course, as 
we have already emphasized, we can, using our values of Co(y max ) from the data, determine 
fi A via Eq. (1291. 



The MH case is somewhat different, particularly with regard to induced singularities 
and instabilities due to dR/dr = C' Q (y) = at z = z max , which in the MH formulation 
now appears in the denominator of the some of the key equations - though not in MHE's 
solution to the null Raychaudhuri equation (122]) . as we pointed out above. R(r) is the MHE 
designation for the angular-diameter (observer-area) distance down our past light cone, for 
which we have used Co(y) in our OC treatment. Once z = z(r) is calculated and the data 
are expressed as functions of r, M(r) must be found by solving the differential equation (see 
MHE, Eq. (14)): 



dr \dR/dr) [dR/dr, 



'dR\ 2 1 1 ~ 
dr +1 ~3 AR 



(39) 



where we have included the A term which MHE set to zero in their treatment. Here we see 
clearly the radial derivative of R in the denominators, which induces a singularity in these 
terms at the maximum of R. Once M(r) is determined, E(r), which in the MH formulation 
gives the local energy per unit mass of the cosmological fluid (which is essentially carries the 
same information as the potential F(y) in the OC formulation), is found via the algebraic 
relationship 



l 

1 + 2E = - 



(#) 2 + i 



(dk\ 2 



M 

-■ (40) 



\dr J 

Obviously, once M(r) is known we can solve Eq. f )29|) with M = M max and C max = R ma x 
as an algebraic equation for A. This complements the MH formulation already in the lit- 
erature by showing how A can be determined in that approach, using the maximum of the 
angular-diameter distance. M(r) and E(r) are two of the three arbitrary functions which 
must be determined in order to obtain the solution in the MH scheme - the "the bang time" 
t#( r ) is the other. In our OC framework, we, instead, focus on the equivalent problem of 

19 



using the data to exactly determine the metric variables as functions of w and y. 



The singularities in Eqs. ( )39|) and ( 1401) . and in other equations of the MH formulation 
of LTB cosmologies, present a challenge to the integration scheme. In their numerical inte- 
gration of the equations Hellaby and Lu [l^] and Hellaby and McClure [15 1 overcome these 
difficulties by carefully matching the solutions on either side of this boundary. As men- 
tioned above, this problem does not arise in the OC formulation. There, as we have seen, 
the equivalent C' (y) = at y max never appears in the denominator. This difference also 
significantly differentiates the two stability analyses as well. In the next section, we present 
a brief analysis of the stability of the null Raychaudhuri equation in its OC formulation. 
As we recall, solution of this equation is the first step in the OC integration process. This 
compares favorably with that in the MH formulation, where again singularities are involved 



15| 



Though there are some good reasons to follow the Hellaby and Lu, and the Hellaby 
and McClure, preference for avoiding fitting the data (e.g. angular-diameter distances and 
galaxy number counts as functions of redshift z) with smooth data functions, there are 
equally good reasons for doing so, as we have briefly indicated earlier. Of course, as in the 
MH approach, we need first to purify the actual data by removing errors due to selection 
effects, absorption, proper motions, image distortions, etc. Once that is done, we find the 
best functional fit to the corrected data using a function which has the proper overall charac- 
teristics that parameter (e.g. angular-diameter distance) must manifest, in terms of limiting 
values (e.g. at z = 0), maxima and minima, inflection points, etc. These are determined 
from the basic theory underlying the definition of the parameter, and the cosmological 
model itself, and severely limit the type of functions we can use. Furthermore, using smooth 
data functions fit to the observed data minimizes problems stemming from redshift-binning 
of data in the fully discrete approach. Finally, as we have pointed out above, having at least 
a formally analytic solution - particularly in simple cases - can illuminate and guide our 
numerical work in more complicated ones, and provide a way of validating it. Now we turn 
to examine the stability of the null Raychaudhuri equation. 
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VII. STABILITY OF THE NULL RAYCHAUDHURI EQUATION 



As uncertainties and errors in the data are introduced into the OC field equations, how- 
are the solutions affected? Will they experience any "blow-ups" from small or modest errors 
or uncertainties? 

We have already seen in the last section that we do not expect any - at least not from 
singularities in the equations themselves. But we now, as an example, look more carefully 
at the first step of the integration procedure, the solution of the null Raychaudhuri equation 
[see Eq. (|22|) above], to see if it is stable against blow-ups induced by small errors or 
uncertainties. We recall that its solution enables us to find z = z(y), which will enable us 
to write other cosmological functions in terms of y instead of z. 

For our purposes in this section it is convenient to write the first integral of the Ray- 
chaudhuri equation, Eq. (1331 . as 



<f>{z,C {z),M {z)) 



A + z) 2 C, z (z) 



A 



l-lf z (l + z)C (z)M (z)dz 



(41) 



where, = dy/dz. The change in due to uncertainties and errors in the data is given by 



A0 



oC oMq 



and, it follows from Eq. fj4T|) that 



(42) 



where, 



dC 



dM 



A I(z)C (z)M (z) 
A I(z)Co(z)M (z) 



C , z (z) 
MqJz) 



Hz) 



A + z) 



2AI 



l-U*(l + z)C (z)M (z)dz 



(43) 
(44) 

(45) 



It is clear from the above that the changes in will be small as long as 5Cq and 5M are 
small. It follows from Eq. (14^1) that, at C max = C(wq, z max ), d(p/dM vanishes. Hence, 



[A(f))z max — AoI max C max M max 5CQ 



(46) 
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where, I max = I(z max ) and M max = M§{z max ). Moreover, it is also clear here that there are 
no singularities that could cause divergences. 

VIII. DETERMINING A UNIVERSE MODEL GIVEN DATA FUNCTIONS 



This section discusses the implementation of the algorithm we presented in [TJJ for using 
the spherically symmetric OC field equations with data to determine the metric. [In Ap- 
pendix |A] we briefly outline the OC integration when we are using redshift-drift data, rather 
than number counts, to gives mass-energy density] We assume that the discrete observa- 
tional data with its gaps and uncertainties are smoothed it into "data functions". For the 
sake of illustration, here we use a so-called ideal data set, that is data functions resulting 
from a smoothing process modeling the Universe as being exactly flat FLRW with A = 0. 

The ideal FLRW data, for the flat case with A = 0, have the form (Stoeger, et al [6]) 

C (z) = 2H^(1 + z)- 2 {z + 1 - (z + 1) 1/2 } (47) 

and 

M (z)=3H (l + z) 1 / 2 , (48) 

where H is the Hubble parameter measured at w = wo and y = 0. 

Solving the null Raychaudhuri equation with this data (see Stoeger aZ [6j) yields the 
following relation between redshift and the null coordinate y 

(1 + *)= v» . ( 49 ) 

a^^ = -. (50) 

From Eq. ([2} and the substitution of Eq. (jl9]l into Eqs. §E§ and (H7|) . we find that on 
our past light cone 



A(w ,y) = A (l - ay) 2 , 
B(w ,y) = A (l - ay) 2 , 

C(wo,y) = Aoy(l-ay) 2 . (51) 



respectively. 
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Application of the algorithm given in ll| to data gives the time derivatives of the metric 
functions A(w, y), B(w, y) and C(w, y) of all orders evaluated on our past light cone w = wq. 
For the ideal FLRW data given by Eqs. (T4"T|) and (T4"g|) the first few terms are 



Mv) 



Mv) 



2a 2 y 



+ C 1 )A (y), 



I - ay 

A&y) + Aya\C x - 2a) 



+ {C 2 



C\ 



2 1 - ay 

B (y) = 2A a(l-ay), 

B {y) = 2A a[(C 1 - a) - ay{C x - 2a)], 

C (y) = 2A ya(l - ay) 

C (y) = 2A ya[(C 1 - a) - ay{C x - 2a)\ 



My), 



(52) 



where, the constants C n are integration constants to be determined as explained in the 
sequel. 



According to the procedure developed in 



ll| we have that A(w, y), B(w, y) and C(w, y) 



are written as a Taylor series in w. Hence, in particular, 



A(w, y) = A (y) + A (y)(w - w ) + -A (y)(w - w ) 2 + ... 
On our world line (y = 0), the above series reads 



(53) 



A{w, 0) =A [1 + C x {w - w ) + l -C 2 [w - w ) 2 + ...] 



(54) 



As explained in Sec. [TT] above, the freedom we have to choose w as any time parameter 
we like along our world line C, is effected by choosing A(w,0). Hence, we clearly see from 
Eq. (|54"|) that the constants C n appearing in the integration procedure are not completely 
arbitrary or independent. We are free to choose these constants as long as they satisfy the 
constraint of being the coefficients in a Taylor expansion of A(w, 0). Next we illustrate this 
important conceptual issue of the observational coordinates formalism by presenting two 
possible choices of the time coordinate w on our world line and their effects on the constants 
C n . 

(i) Let us choose w to measure proper time along our world line C. That choice [see 
Appendix [B] for details] leads to 
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A(w,0) = 1 



(55) 



and from Eq. (1541) we find that A = 1, and all the constants C n vanish. 

(ii) Consider the usual way in which the flat FLRW (A = 0) metric is written in obser- 

n n 

vational coordinates 6fl, I7J: 



Hence, Eq. Q53} gives C\ = 2a, C 2 = 2a 2 , C n = for n > 3 and A ^ 0. We note that 
the series is truncated at the third term. Therefore, from the expressions fl5~T]) and (152]) we 
can reconstruct the solution corresponding to the specific choice of time coordinate on our 
world line given by Eq. (15"T|) . 

IX. CONCLUSION 

Here we have focused on key theoretical and mathematical aspects of determining the 
large-scale metric of the universe from cosmologically relevant astronomical observations, 
comparing and contrasting the traditional MH 3+1 LTB and the observational cosmology 
(OC) approaches to formulating and integrating the inhomogeneous spherically symmetric 
field equations with a cosmological constant, A ^ 0. Such approaches are crucial, since, 
despite the attractiveness and success of almost-FLRW models, there are significant uncer- 
tainties about whether or not they provide the unique best-fit description of the large-scale 
structure of the cosmos. These uncertainties can really only be resolved by non-perturbative 
treatments which go beyond FLRW. In pursuing such paths we may find that we can con- 
firm that, yes, the universe is indeed almost-FLRW on the largest scales, and that the 
Copernican principle holds. In this case, we would also be able to determine more precisely 
the smallest length-scale on which the universe is almost-FLRW (it is clear from its small 
and intermediate scale lumpiness that it is not almost-FLRW on those scales). However, 
a second possible outcome of such studies is that the universe is not almost-FLRW on the 




(56) 



It follows from Eqs. (136"]) and (150]) that 



A(w, 0) = A (aw) 2 



(57) 
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largest scales, or was not or will not always be so. These approaches, beginning with spher- 
ically symmetric inhomogeneous models, would then provide a more secure best-fit model 
to the observations, and a measure of how much deviation there is from the cosmological 
principle. More importantly they would, as we have discussed in this paper, confirm whether 
or not there is indeed a significant amount of vacuum energy, or dark energy, in the universe 
- or whether instead the apparent acceleration of cosmic expansion is due to large-scale 
inhomogeneities. We may end up confirming that the universe is almost-FLRW with 73% 
dark energy, as the concordance cosmic model says, or we may find that it is really close 
to an inhomogeneous LTB with or without dark energy. Either way, these programs would 
significantly contribute to our cosmological knowledge and understanding, and our level of 
confidence in it. 

In the course of the paper we have summarized the essentials of both the OC and MH 
approaches - the first employing observational coordinates suited to the past-light-cone 
structure of our observations made here and now on our world line, and the second relying 
on a 3 + 1 foliation of space-like surfaces orthogonal to the cosmic time coordinate. In 
particular we sketched the essential steps in integrating the field equations with data in 
each of the two formulations, high-lighting the main similarities and the differences. In the 
OC case we also show - in Appendix [A] - how redshift drift can be used to give us the 
mass-energy density of the universe as a function of redshift - once we are able to obtain 
such data - instead of the more problematic galaxy counts (this obviously can also be done 
in the MH formulation). Although both formulations and integration procedures promise 
to lead to the same results, the OC procedure manifests certain advantages - particularly in 
the avoidance of coordinate singularities at the maximum of the angular-diameter distance, 
and in the stability of the solutions obtained. The OC procedure is, as already emphasized, 
geared towards determining the metric variables from observations, whereas MH procedure 
is directed towards determining the three free functions that appear in that formulation. 
These outcomes are equivalent, of course, and each perspective has advantages. But the 
OC approach enables one to link the metric more directly to observations. Finally, the 
possibility of performing some analytic integrations, at least in simple cases, by first fitting 
the corrected data to smooth data functions seems to have definite advantages to simply 
numerically integrating the difference equations directly with binned data. Smoothed data 
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functions enable us to include properties that data must have within the model. 



This, together with the OC procedure itself, is illustrated in Sec. IVIIII with the very 
simple abbreviated example of deriving the FLRW metric from FLRW data. We have 
treated this before [?| by obtaining the solution on our past light cone. Here we completed 
that solution by moving it off the light cone, as specified in Sec. [V] We should note that, 
although FLRW data functions can be derived from the FLRW metric, it is non-trivial - 
and must be demonstrated - that FLRW data functions uniquely imply the FLRW metric. 
We suppose that we have ideal FLRW data, including knowing that our universe is flat and 
that z max = 1.25. We then carry out our integration scheme for that data - showing that 
it leads to a flat FLRW metric with A = 0. This demonstrates the advantage of having 
data functions which can be used to construct analytic solutions to the field equations. 
No singularity problems were encountered in carrying out this integration. We could have 
included uncertainties or errors in the data functions, of course, and then determine how 
they propagate into the solution, and whether or not they are large enough to affect its 
stability. This will be pursued later. 

In a future paper we shall be investigate how cosmic microwave background (CMB) and 
CMB anisotropy observations fit into this OC integration scheme, and to what extent they 
give us information independent from other sources. 



Appendix A: 

First, we briefly derive a relationship which enables us to find z = z(y) using redshift-drift 
data - without having to solve the null Raychaudhuri equation, Eq. ( !22l . 

It is clear that on the radial null geodesies r = r(w,y) and z = z(w,y). Now, since on 
our past light cone w = Wq - in OC coordinates w labels past light cones - we must have 
that 

dr Adw — Bdy (Al) 

dz zdw + ^dy ' 

where, z = dz/dw is the redshift-drift. It follows that dividing the numerators and denomi- 
nators on the right-hand side of Eq. (I A II) by dw gives dr/ dz = A/z. If, instead, we divide the 
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numerators and denominators on the right-hand side by dy we obtain dr/dz = —B(dy/dz). 
Equating these two results gives 



B(w Q ,y) Idz 



(A2) 



A(w ,y) zdy 

[See [12] for a detailed derivation on this result.] Now, using the freedom to choose the radial 
null coordinate y by setting B(wo,y) = A(wo,y) [see Sec. [II], Eq. (1A2j) is further simplified 
on our past light cone. 

Since z(z) is given from data, solving Eq. (1A2I) gives z = z(wo,y). That is the same 
information that would be obtained from the null Raychaudhuri equation f l2"2"|) . It is quite 
clear that, whereas to solve the latter on our past light cone one needs the mass-energy 
density fi(wo,z), what we need to solve Eq. (1A2|) is a different piece of information - the 
redshift-drift z(w ,z). 



Secondly, Araiijo and Stoeger [12] have shown that using Eq. flA2p and its solution, we 
can write the mass-energy density in terms of the redshift z as 



2zp(z) d 
M(z) dz 



1 + Zq(z) 



-4u)n(z). 



(A3) 



Eq. ( 1A3I) shows that the mass-energy density fi mo (z) can be completely determined in 
terms of the redshift z, the redshift-drift z(z) and observer area distance C(wq,z) data 
on our past light cone. We observe that the ^ mo {z) dependency on the observer area dis- 
tance C(wq, z), leads to its dependency on /ia, that must be determined by data giving the 
maximum of the observer area distance, C (wo, 2 max ), and the redshift z max at which that 
occurs. 

Now, from Eqs. ©, (jAJ) and (USD we find that 



u Q (z) 



1 



1 + 



z dC (z) 



2Cn(z) - Z 



dCp(z) 
dz 



(A4) 



2Cl(z) \ Al{z) dz 
Clearly, dC (z)/dz can be determined from the r (z) = C(w , z) data, through fitting. We 
determine Cq{z) by solving Eq. ( 124"1) for Co(y) on w = w [see Araiijo and Stoeger for 
details] and then use Eq. ( 1A2j) and its solution to write the result in terms of z. Taking into 
account the appropriate boundary condition (central condition), Cq(z) is given by: 
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C (z) 



2C (z)Jo Zq(z) 



Am 



z 2 (z) 



.^dC (z) 



dz 



2C (z)dA (z) | dC (z) 
A (z) dz dz 



A 2 (z)C 2 (z)fi A \dz. 



(AS) 



Since, z(wq,z) is very small, it follows very clearly from Eqs. (1A3|) . ( ]A4j) and ( 1A5j) that 
a useful approximation to fi mo (z) is 



// mo (z) S -Acuo(z). (A6) 

First, we note that having determined n mo (y) from data, we can easily find M(y) by 
Eq. (1521) . once the potential is known. But, F(y) is completely determined by data 
through Eqs. ( 1271) and ( 1A5I) . except for the unknown value of the constant yUA, which we 
carry along for the time being. Then we evaluate the mass function M(y) at y ma x and plug 
the result into Eq. (|29|) . which with the C max data becomes an algebraic equation for fi^. 
With this determination of fi\, we know Co{y) completely, and can now determine F(y) 
more precisely from Eq. (|27l) . Furthermore, we observe from Eq. ( fl~9l) that the quantity 
Uo(y) is also completely determined at this stage. We are now ready to following our solution 
off out past light cone w = wo for all w applying the procedure described in Sec. [V] from 
Eq. (135!) onwards. 



Appendix B: 

Ellis et al. jlj] have shown that if the coordinate w is normalized by the condition that it 
measures proper time along C, and r denotes proper time along general galactic world lines, 
then the redshift observed at C for the light coming from a galactic world line is given by 

l + z = - (Bl) 

with the ratio dw/dr evaluated along the galactic world line. 

Since y is a comoving coordinate, and dr 2 = —ds 2 , then Eqs. ([1]) and (1B1|) lead to 

dw 1 
Z ~ dr ~ A(w,y) 
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where the normalization at y = (our world line) gives 



>Lo~^°) _i (B3) 
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